Bidirectional sorting of flocking particles in the presence of asymmetric barriers 
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We numerically demonstrate bidirectional sorting of flocking particles interacting with an array of 
asymmetric barriers. Each particle aligns with the average swimming direction of its neighbors ac- 
cording to the Vicsek model and experiences additional steric interactions as well as repulsion from 
the fixed barriers. We show that particles preferentially localize to one side of the barrier array over 
time, and that the direction of this rectification can be reversed by adjusting the particle-particle ex- 
clusion radius or the noise term in the equations of motion. These results provide a conceptual basis 
for isolation and sorting of single- and multi-cellular organisms which move collectively according 
to fiocking-type interaction rules. 

PACS numbers: 87.10.-e,05.65.+b,87.17. Jj,05.40.-a 



I. INTRODUCTION 

The ensemble dynamics of self-driven particles can dif- 
fer significantly from those of Brownian random walkersl . 
For example, in experiments on microfabricated habitats 
connected by funnel-shaped channels, self-propelled E. 
coll bacteria preferentially migrated to the chamber to- 
wards which the funnels pointed^i^ even though Brow- 
nian particles would have remained equally distributed 
in both chambers. A simple simulation model showed 
that the rectification arises due to the modification of 
the run-and-tumble swimming dynamics of the bacteria 
by the walls of the microenvironmenf*. When a run- 
ning bacterium encounters a wall, it does not reflect away 
from the wall or immediately tumble, but swims in the 
direction of the wall while preserving as much as pos- 
sible its prior direction of motion. Refsj^ and^ found 
rectification under this interaction rule for independent 
swimmers that did not interact with each other. The 
rectification in the bacteria system resembles a ratchet 
effect in which a net dc motion occurs in the absence of 
a dc drive due to the application of an external ac drive 
or flashing substrate^. For self-driven particles, however, 
no external driving is necessary. In addition to demon- 
strations of directed bacterial motion achieved through a 
ratchet mechanisroi, it has also been shown that baths of 
swimming bacteria can induce directed rotational motion 
of asymmetric flywheels^"^''. 

Interactions between self-propelled particles can lead 
to distinctive dynamical behaviors that are more com- 
plex than those of independently moving particles. Sim- 
ple models such as that of Vicsek et alJ^^^ capture 
many features of the dynamics of species with strongly 
collective motion, in which individuals preferentially 
align with their neighbors and form moving groupsi^ii^. 
These models qualitatively reproduce the motion of 
both macro-scale groups, such as fish schools and bird 
flock a^^i^^ , and micro-scale groups, such as bacterial 
swarms and cancerous tumors^^. The original Vicsek 
model includes only a term for preferential velocity align- 
ment with all neighbors within a fixed flocking radius. 



yet it exhibits a phase transition to unidirectional mo- 
tion as a function of particle density and noise ampli- 
tude. Although numerous modifications of the Vicsek 
model have been proposed, such as the addition steric 
interactions—- and/or cohesionl2r=, only a very limited 
amount of work has been done on the interaction of flock- 
ing particles with walls or barriers. Walls can impose a 
directional symmetry breaking^'*, induce the formation 
of a vortex stat o^^'^^" — or laningS^, or simply serve as 
aggregation focal points^^; walls have also been used for 
understanding finite size effects^Zr— , such as the relation- 
ship between the collective dynamics of fish in a tank and 
those of fish in the open ocean. 

In this work we simulate a modified version of the 
Vicsek fiocking algorithm that includes both steric re- 
pulsion between particles and confinement within a two- 
dimensional microenvironment with strategically placed 
gates similar to those of Refi^. Here we consider strictly 
repulsive particle- wall interactions, so that the particles 
do not follow the walls when swimming independently. 
As the particle density increases, we find rectification ef- 
fects once the density is high enough to permit collective 
motion to occur. By varying the interparticle exclusion 
radius, the flocking radius, or the noise, we can reverse 
the direction of the rectification. This result has implica- 
tions for the potential sorting of self-propelled particles 
that move according to these types of interaction rules. 

Simulation- We consider a two-dimensional L x L sys- 
tem of TV self-driven particles at number density po = 
N/L/^ with fixed, repulsive boundaries on all sides. The 
overdamped equation of motion for a single particle i is 
dxi = Wi(t)dt with 



v,(i) = f:,(i) + C(t) + f^(t) 



(1) 



All quantities are rescaled to dimensionless units. The 
time step dt — 0.002 and we take L — 66. The velocity 
consensus force f^^, also called the alignment force, is 
determined by the velocities of all M particles, including 
particle i, within a flocking radius r/ of particle i: 



iUt) = Af (cos($;,(t))x + sin($;,(t))y) 



(2) 
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Figure 1: Simulation images. Lines: barriers and walls; dots; 
particle positions, a) Initial state of sample with r/ = 1.0 
and = 0.07 at po = 0.4. b) The same sample after 7 x 10^ 
simulation time steps shows rectification of particles into the 
top chamber. 



with 



^Icit) = arctan' 




(3) 



Here Af = 2.0 and ^ is a random variable uniformly 
distributed on the interval [—77/2,77/2]. Both the steric 
particle-particle interactions f* and the particle-barrier 
interactions are given by the stiff spring repulsions: 



Ar(2re - rij)Q(2re - rij)rij and f^(i) = 

'rik)Q{re + rg-rik)rtk, where Ar = 200, 
10, Tij = \ri{t)-rj{t)\, and r^j = [ri{t)-rj{t)]/rij. 
Here Tg is the particle exclusion radius, rg — 0.05 is the 
barrier exclusion radius, and there are Ng = 16 barriers 
composed of the four confining walls plus 12 V-shaped 
gates. is the vector from the nearest point on barrier 
k to particle i, nk = \r^k\, and r^k = T^ik/nk- The length 
of each side of the V gates is Lb = 4.9 and the angle 
each V arm makes with the y axis is 30°. The spacing 
between the bases of the V's is Is = 5.5 and the spac- 
ing between the tips of adjacent V's is — 0.6. The 12 
gates bisect the system into top and bottom chambers, 
with the aperture of each funnel shape pointing toward 
the top chamber. We initialize the system by distribut- 
ing the particles at random throughout the sample. The 
equations of motion are then integrated for 3 x 10^ sim- 
ulation time steps. In the absence of particle-particle 
interactions, the purely repulsive wall interactions pro- 
duce no rectification of the particles, in agreement with 
the results of Ref.'^ 

Results- In Fig. [ija), we show an image of the sim- 
ulation geometry in the randomly initialized state for a 
system with rj = 1.0, = 0.07, and po = 0.4. After 
a sufficient amount of time elapses, the particles concen- 
trate in one of the two chambers, reaching a steady state 
value of ptop, the density in the top chamber. In Fig.[ljb), 
after 7 x 10^ simulation time steps the particle density is 
clearly higher in the top chamber. 

We find that we can vary whether the rectification 
moves the particles into the top {ptop > Po) or bottom 
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Figure 2: (a-c) ptop, the density in the top chamber, after 
3 x 10" simulation time steps for a sample with initial density 
Po = 0.4, indicated by the dashed line, (a) ptop vs for 
»7 — 1.5 and r/ = 1.0. The rectification reverses at Ve = 0.12 
and drops to zero for > 0.3. b) ptop vs i] for = 0.12 and 
r/ — 1.0. The rectification reverses at 77 ~ 1.0 and drops to 
zero for rj > tt. c) ptop vs r/ for — 0.12 and 77 — 1.5. For 
Vf < Te only steric particle interactions occur and rectification 
is negligible. There is a rectification reversal at r/ ~ 1.2, 
and for large r/ when all the particles tend to align into a 
giant flock, the particles accumulate in the top chamber, d) 
Rectification phase diagram for vs rj. Lower contours (red) 
indicate rectification into the top chamber and upper contours 
(blue) indicate rectification to the bottom chamber. 



(ptop < Po) chamber by altering ry, rg, or 77, as shown 
in Fig. [21 a-c) where we plot ptop after 3 x 10^ simula- 
tion time steps. For small values of re and 77, particles 
are rectified into the top chamber, but a rectification re- 
versal occurs at = 0.12 and 77 ^ 1.0 in Fig. [2ja) and 
(b). respectively. There is a saturation into a nonrecti- 
fying state for re > 0.3 in Fig. [5fa); this corresponds to 
2re > lo and occurs when the particle diameter becomes 
larger than the aperture between adjacent gates, so that 
particles can no longer pass between the upper and lower 
chambers. In Fig. [H^b) rectification vanishes for 77 > tt 
when the alignment force between neighboring particles, 
and thus the tendency of particles to form fiocks, is al- 
most completely destroyed by noise. We plot a rectifica- 
tion phase diagram as a function of re and 77 in Fig.[2I^d). 
showing that rectification into the upper chamber occurs 
for small values of re and 77, while reversed rectification 
into the lower chamber appears for larger re and small 77. 

We next consider the fiocking radius r/. Fig. ^c) in- 
dicates that no rectification occurs when rj < re- In 
this limit, the particles interact only sterically and have 
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Figure 3: Illustration of rectification into the top chamber 
for low noise 77 = 1.5 and small exclusion radius — 0.05 at 
r/ = 1.0. A 15 X 15 section of the sample is shown. Dots: par- 
ticle positions; light lines: particle trajectories; heavy lines: 
barriers. A flock incident on the gates from the bottom cham- 
ber (a) condenses and elongates in order to file through the 
aperture between gates (b). Flocks incident on the gates from 
the top chamber have a much lower probability of passing 
through the aperture and cannot be tunneled into a similar 
oblong shape. 



no flocking interaction, and the repulsive barrier walls 
produce no rectification in the absence of fiocking. For 
Te < rf < 1.2, we find a reversed rectification into the 
lower chamber, while for all r^^ > 1.2, the particles rec- 
tify into the top chamber. For > 2.0, the value of ptop 
saturates at ptop — 0.8 — 2po, indicating that nearly all 
of the particles are located in the top chamber. 

The rectification reversal occurs due to a change in the 
nature of the microscopic interaction between the flocks 
and the funnel channels. For example, as the exclusion 
radius increases, the particles are less able to form 
tight and cohesive flocks. At low values of Vg, particles 
are rectified into the top chamber when flocks, incident on 
the gates from the bottom, rearrange into oblong shapes 
and pass efficiently through the funnel, as illustrated in 
Fig. [31 For higher values of r^, the steric interparticle 
repulsion prevents the flocks from condensing and makes 
it impossible for more than one particle at a time to pass 
through the funnel aperture. As a result, the particles 
clog inside the funnel rather than passing through, as il- 
lustrated in Fig. HJa). The flock reverses direction due 
to the repulsion from the barrier walls, and at most one 
or two particles occasionally manage to escape the flock 
and enter the top chamber, as shown in Fig. IHb). In 
contrast, a flock that approaches the gates from the up- 
per chamber is fragmented by the gates into two smaller 
flocks; when this occurs, particles that are directly inci- 
dent on the aperture between gates can escape from both 
flocks and pass in a single file into the lower chamber, as 
illustrated in Fig. IHc,d,e). Since the average number of 
particles escaping the flock and crossing the barrier is 
larger when the flocks are approaching from above than 
when they are approaching from below, a net rectification 
into the lower chamber occurs over time. We note that 
the reversed rectification into the lower chamber (Fig. 0]) 
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Figure 4: Illustration of rectification into the lower chamber 
for low noise rj — 1.5 and large exclusion radius re = 0.2 at 
r/ = 1.0. A 20 X 20 section of the sample is shown. Dots: par- 
ticle positions; light lines: particle trajectories; heavy lines: 
barriers. (a,b) Flocks incident on gates from the bottom 
chamber cannot fit through the aperture for this value of Ve; 
the flock jams inside the funnel while a single particle (high- 
lighted in red) escapes from the flock and enters the upper 
chamber. The remainder of the flock returns to the lower 
chamber. (c,d,e) Flocks incident on gates from the top cham- 
ber are fragmented and a small group of particles (highlighted 
in red) can escape from the flock and enter the lower chamber. 
The flock fragmentation process occurs with greater frequency 
as the flocks become less cohesive due to either higher rj or 
higher ratios r^/rf. 



is a much slower process than the forward rectification 
into the higher chamber (Fig. [3]), although we are able 
to reach a steady state within our simulation time for ei- 
ther process. In spite of this, we find that the maximum 
possible amount of rectification that can be achieved in 
steady state as the parameters are varied is the same for 
both directions of rectification, as shown in Fig. [2{a). 

Reversed rectification into the lower chamber also oc- 
curs whenever the flocks become fragile or prone to break- 
age. This occurs both when the noise parameter 77 is in- 
creased and when the flocking radius is reduced. Un- 
der these conditions, the flocks are not cohesive enough 
to flow as a unit through the funnel aperture in the man- 
ner illustrated in Fig. [31 at the same time, the probability 
that a flock will fragment and lose some of its members 
to the lower chamber when approaching the gates from 
above, as in Fig. Hl^c-e), is increased. In Fig. [l^a), we 
plot the average flock size Nc as a function oi rf, r^, and 
7/ for the systems in Fig. [2l We separate the particles 
into clusters iteratively by identifying particles that are 
within the flocking radius r/ of each other; Nc is then 
the average number of particles per cluster. The value 
of Nc is higher in regimes where the particles are recti- 
fled to the top of the container, and lower in the reversed 
rectiflcation regime. 

One of the unique aspects of the rectification behavior 
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described here is that, unhke previous rectification phe- 
nomena reported for self-driven particles^, it occurs only 
when the initial particle density po is high enough for 
flock formation to occur. In the limit of low poi when 
the particles are moving independently and not able to 
form flocks, individual particles simply reflect off the bar- 
riers in a manner similar to inertial particles. This type 
of barrier interaction has been shown to produce no rec- 
tification in noninteracting particle limit^, and as indi- 
cated in Fig.[6^a) we find no rectification at low densities 
Po < 0.01. As Po increases, both rectification and a recti- 
fication reversal emerge, and the amount of rectification 
saturates for po > 0.1, as shown in Fig. ^h). We note 
that since po represents number density, rather than sur- 
face area covered, it is possible to have po > 1. 



Figure 5: Mean flock size Nc, in number of particles, vs r/ 
(blue circles), (red squares), and rj (green diamonds). Error 
bars indicate standard deviation. 
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Conclusion- We have implemented a simple model of 
flocking particles in the presence of flxed, repulsive bar- 
riers, and flnd that such particles will concentrate on one 
side of a set of asymmetric V-shaped gates. The direc- 
tion of the rectiflcation can be reversed by modulating 
any of three parameters: the flocking radius r/, the ex- 
clusion radius r^, or the noise parameter rj. The exis- 
tence of the rectification and its direction are determined 
by the ability of the particles to form flocks and the ro- 
bustness of the flocks against breakage; in the low den- 
sity limit, when no flocks appear, we find no rectification 
due to the purely repulsive interactions of the particles 
with the barrier walls. Thus, the rectification we observe 
arises strictly due to collective effects. The bi-directional 
rectification behavior we describe could be used to sort 
particles which tend to concentrate on different sides of 
the barrier?2.. This effect is similar to the sorting phe- 
nomenon observed by Mahmud et al. for cancer cells^i. 
We expect sorting devices based on these principles to 
have broad potential applications with regard to both 
biomedical diagnostics and therapeutics. 



Figure 6: Dependence of rectiflcation on initial particle den- 
sity Po for a system with r/ — 1.0 and iq — 1.1. (a) ptop/ po 
after 3 x 10^ simulation time steps vs re for difi'erent values 
of po. From blue to red, po = 0.004 (o), 0.01 (□), 0.03 (x), 
0.05 (A), 0.08 (v), 0.1 (O), 0.12 (+). A rectiflcation reversal 
emerges as po increases, b) ptop/ po after 3 x 10® time steps vs 
po for re = 0.05 (upper red curve) and Ve ~ 0.23 (lower blue 
curve) . 



This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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Figure 7: Supplemental Figure. Demonstration of two-species 
sorting. Simulation of system containing 871 "A" particles 
with re = 0.22 (green) and 871 "B" particles with = 0.055 
(pink), rf = 1.0 and rj = 1.1 in both cases. Particles of dif- 
ferent species repel via steric repulsion but only experience 
alignment forces with particles of the same species. Simula- 
tion shown (a) at time t — and (b) after 4 x 10^ simulation 
time steps when many of the "A" particles have rectified into 
the bottom chamber and most of the "B" particles have rec- 
tified into the top chamber. The bidisperse system requires 
longer times to reach a steady state compared to the monodis- 
perse system. 
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